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Proximate and ultimate factors that promote aggregated 
breeding in the Western Sandpiper 


Matthew Johnson’, Jeffrey R Walters 


(Department of Biological Sciences, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061, USA) 


Abstract: We report that Western Sandpipers (Calidris mauri) on Alaska’s Yukon-Kuskokwim River Delta 
exhibited aggregated breeding behavior at a relatively small spatial scale. Prior to clutch initiation, males performing song 
flight displays on a 36 ha plot were aggregated as were subsequent initial nesting attempts on the plot. We tested three 
hypotheses commonly invoked to explain aggregated breeding in territorial species (social mate choice, predation, and 
material resources hypotheses), and found support for the material resources hypothesis, as dispersed individuals were 
more often associated with tundra habitat patches, and aggregated individuals nested more often in undulating-tundra 
habitat patches compared to patch availability. The pattern of habitat occupancy conformed to an ideal despotic 
distribution with aggregated nesting birds in undulating-tundra patches experiencing lower reproductive success. On our 
study plot, older, more aggressive males solicited females more often, and defended larger, more dispersed sites in tundra 
habitat patches, compared to younger, less aggressive males that were aggregated in undulating-tundra habitat patches. 
Breeding aggregations are often concentrated on or near a critical resource. In contrast, Western Sandpiper breeding 
aggregations occur when dominant and/or older individuals exclude younger, subordinate individuals from preferred 
habitat. Although many taxa of non-colonial birds have been reported to aggregate breeding territories, this is the first 
quantitative report of aggregated breeding behavior in a non-colonial monogamous shorebird species prior to hatch. 
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BRIM. EIRP, FKEA ER Berk PE FIRE SS Ta) HE SB EI BB 2, JE HAR TS Ai SASF Le DE 
ERE ROK ISU FALE F, AES AB PE BOS RE S WU RE ASS OR A JA SE UTP, SRE 
I RARE 2 SRB PR SBE DRL A A Bd is DBASE, 1 i mi SI HP) H BU a h ER Ach Ft AA 
MREFA FME HA AS PAE Je EE JE EAP Jd» EER, SPV Sa tah YS AE Se EAR ERE ee) SZ Ha 
FL FRIŠ, TATA SC HE PS CSTE ARE JE A SEP aif SS IF ZB TS AP AE AT SET FE 


ABE ial: VEYA; ANTAL; PR; SE) A; DE BR Ha AB 
FARIS ES: Q958.155;Q959.722 XERRED: A MBAS: 0254-5853-(2011)02-0128-13 






















































































= 


























r 




































































= 
at 
























































[= 












































zN 

















(s3 

































































































































































H 



































Nil 
























































jrs; 



















































































= 












































Breeding aggregations have been observed in (Itzkowitz, 1978), insects (Muller, 1998; Sumpter & 
numerous aquatic and terrestrial animals including fish Broomhead, 2000), reptiles (Stamps, 1988), primates 
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(Treves, 2000), and birds (Stamps, 1988; Muller et al, 
1997). Presumably, animals aggregate when the benefits 
of doing so exceed the costs of social interaction. The 
benefits of joining a breeding aggregation may be either 
resource-based, where (in birds) individuals select 
habitats that offer foraging opportunities or safe nesting 
sites (natural selection), or driven by sexual interactions 
and mate choice (sexual selection hypotheses, Tarof & 
Ratcliffe, 2004; Tarof et al, 2005). Mechanistic and 
evolutionary processes that result in avian breeding 
aggregations have most often been studied among 
colonial nesting species where nesting and is the primary 
activity within the aggregation, and among lekking 
species where only mate choice and reproduction occur 
(reviewed by Bradbury & Gibson, 1983, Brown & 
Brown, 2001). Among colonial nesting species, potential 
costs of nesting in close proximity to conspecifics 
(increased levels of: resource competition, intraspecific 
brood parasitism, nest predation, and parasite transmission) 
appears to be offset by direct benefits including 
increased food finding efficiency (material resources 
hypothesis, Kiester & Slatkin, 1974; information center 
Wittenburger & Hunt, 1985), 
predation due to group defense or the dilution effect 


hypothesis, reduced 
(predation hypothesis, Hamilton, 1971), or increased 
opportunities for extra-pair copulations (hidden lek 
hypothesis, Wagner, 1997). Among lekking species, 
theory on the evolution and maintenance of aggregative 
behavior suggests individuals aggregate in response to: 
the patchy distribution of ecological factors (e.g. display 
sites or predators, material resources and predation 
hypotheses), the patchy distribution of conspecifics 
(males cluster in areas where females tend to be found, 
low-quality males cluster around high-quality males), or 
because females prefer aggregations of males (social 
mate choice hypothesis; Allee, 1951; Darling, 1952). 
Aggregation of multi-purpose breeding territories is 
another common form of spatial structure in animals, and 
occurs when individuals establish territories in close 
proximity to conspecifics which results in clusters of 
territories in apparently homogeneous breeding habitat 
(Tarof & Ratcliffe, 2004). Study of breeding aggregations 
among socially monogamous avian species (>90% 
worlds’ species; Lack, 1968) has revealed that the spatial 
distribution of breeding activities is most often correlated 
with optimal displaying/nesting habitat, areas with high 
food availability, low predation risk, and species specific 
foraging and reproductive strategies (material resources, 


predation, and hidden lek hypotheses; reviewed by 
Hildén(1965), Cody (1985) and Morse (1989); but see 
Tarof & Ratcliffe (2004)). Monogamous avian breeding 
aggregations offer a manageable natural system for 
testing alternative hypotheses explaining the relative 
influence of social and environmental factors on the 
spatial distribution and behavior of individuals. Avian 
breeding behavior is typically observable, aggregations 
can be measured accurately and habitat characteristics 
can be quantified with reasonable effort. Shorebirds 
(Charadriiformes) exhibit some of the most elaborate and 
overt breeding displays of any avian group (Miller, 1985). 
The 24 species of sandpiper in the subfamily 
Calidridinae (Charadriiformes, Scolopacidae) exhibit 
every described mating system, with territoriality and 
mating strategy sometimes varying within a species 
(Oring, 1982). Studies of intra- and inter-specific 
variation in calidridine sandpiper social systems have 
contributed greatly to understanding the evolution of 
social behavior (Pitelka et al, 1974; Miller, 1979; Oring, 
1982; Oring & Lank, 1985). The study of calidridine 
sandpiper spatial structure offers an opportunity to extend 
our understanding of how ecological and social factors 
affect the distribution of breeding animals. 

During six years of study in Western Alaska 
(1999-2001, 2003-2005), we observed apparent 
aggregation of initial nesting attempts by a territorial 
calidridine sandpiper, the Western Sandpiper (Calidris 
mauri, Fig. 1). Western Sandpipers are socially and 
genetically monogamous (i.e., extra-pair paternity is rare) 
and exhibit biparental care of eggs and young (Holmes, 
1971, 1973; Blomqvist et al, 2002). After arriving on the 
breeding grounds, males establish territories on upland 
tundra habitat that are defended against conspecifics, 
from which territorial males advertise for mates using 
species-specific displays (Brown, 1962; Holmes, 1971). 
Territories range from 0.2—0.3 ha, and males often 
defend display territories in close proximity to one 
another. Territorial males regularly engage in chases and 
intense fighting; however, after some clutches are 
completed, it is not uncommon for unpaired males to 
display over nesting habitat in which females have 
previously initiated nests (Holmes, 1971; Lanctot et al, 
2000). Breeding displays are often performed with one or 
more other males simultaneously, but not in a duetting 
fashion (Lanctot et al, 2000). Adult Western Sandpipers 
obtain some food in the vicinity of their nests but more 
commonly forage in separate feeding areas along the 


130 Zoological Research 





© 2003 initial nest attempts 
REZ] WET MEADOW 














TUNDRA-GRASS 
B UNDuLATING-TUNDRA 
























M Yy 


L? 


gf e Uj 


Fig. 1 Habitat map and initial nesting attempts of Western 
Sandpipers (2003) at the Yukon Delta National 
Wildlife Field Station, 
Yukon-Kuskokwim River Delta, Alaska (2003). 
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Habitat in white areas was not quantified. 


shores of lakes, rivers, and sloughs, and in low-lying 
marshes (Holmes, 1971; Lanctot et al, 2000). 

In our study, we provide statistical support that 
Western Sandpipers exhibit aggregated breeding 
behavior. We then evaluate three hypotheses commonly 
invoked to explain aggregated breeding in territorial 
animals: (1) The social mate choice hypothesis, which 
predicts that aggregated males have improved success at 
attracting social partners compared to solitary individuals; 
(2) The predation hypothesis, which predicts that 
aggregated breeding reduces predation via antipredator 
strategies; and (3) The material resources hypothesis, 
which predicts that individuals aggregate in response to 
patchily distributed resources (vegetation and/or food). 
We tested the social mate choice hypothesis by 
comparing clutch initiation dates and female age between 
aggregated and dispersed nesting birds, the predation 
hypothesis by comparing nest success between 
aggregated and dispersed nesting birds, and the material 
resources hypothesis by comparing vegetation associated 
with aggregated and dispersed nesting birds. 
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1 Materials and Methods 


1.1 Study site 

We studied Western Sandpiper breeding behavior 
at the Yukon Delta National Wildlife 
Kanaryarmiut Field Station, Yukon-Kuskokwim River 
Delta, Alaska (61°22'N, 165°07'W). 
surrounding Kanaryarmiut Field Station is a complex of 


Refuge’s 
Vegetation 


lakes, sloughs, wetlands, and upland tundra corresponding 
to the “lowland moist low scrub” community described 
by Jorgenson & Ely (2001). Upland tundra vegetation is 
dominated by lichens, Sphagnum spp., Betula nana, Salix 
fuscescens, Ledum decumbens, Empetrum nigrum, Rubus 
chamaemorus, Eriophorum spp., and Poa eminens, and 
wetland and lakeshore margin vegetation is predominately 
characterized by sedges (Carex mackenziei, C. rariflora), 
(Poa Calamagrostis spp., 
Eriophorum spp.; Ely & Raveling, 1984). The upland 


and grasses eminens, 
tundra vegetation community is a mosaic of patches that 
can be classified into three basic types (Johnson & 
McCaffery, 2004, Johnson et al, 2009). First is upland 
tundra that contains graminoid species (hereafter tundra- 
grass). The second type is upland tundra that lacks a 
graminoid component (hereafter tundra). The third is 
undulating-tundra, characterized by greater vertical relief 
that results in retention of water for longer periods after 
spring thaw and a sedge and grass component. 
1.2 Field surveys 

The breeding population of Western Sandpipers at 
Kanaryarmiut Field Station was continuously monitored 
for eight years (1998-2005) as part of a long-term 
demographic study; however, effort varied across years. 
In 1998 and 2002, demographic data were primarily 
gathered on a 16 ha plot (1998: 43 nests monitored, 53 
adults banded; 2002: 58 nests monitored, 28 adults 
banded). During the other six years of study at this site 
(1999-2001 and 2003-2005), data were collected on an 
additional 20 ha surrounding the original plot “36 ha 





total; mean + SD number of nests monitored per yr = 
(113 + 25), mean + SD number of adults banded per yr = 
(89 + 16)”. Further, individual behavioral data were 
recorded on the 36 ha plot during the last three years of 
study (2003-2005). We only used data gathered during 
those years when the larger study plot was surveyed 











because of temporal variation in effort and to maintain 
consistency in analyses across years. Two to four 
observers surveyed the 36 ha study plot daily from early 
May through late July for banded birds, nests, and broods. 
Adults and chicks were marked with a U.S. Geological 
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Survey identification band as well as with unique UV- 
stable color band combinations at the nest. The location 
and behavior of banded birds was recorded daily 
(2003-2005). The locations of nests and banded 
individuals were mapped using ArcMAP (ESRI GIS and 
mapping software) and nests were monitored through 
hatch, predation, or abandonment. We produced 
vegetation maps in the field by manually delineating 
vegetation patches (tundra, tundra-grass, undulating- 
tundra) on aerial photographs. Polygons were then 
digitized and georeferenced using ArcMAP. 
1.3 Measures of aggregation 

Any classification of random, aggregated, and even 
distributions in nature comprises arbitrary distinction 
made along a continuum of three spatial patterns. 
However, the homogeneous Poisson process is the only 
point process that is the equivalent of complete spatial 
randomness (Poole, 1974, Pielou, 1977). Testing for 
complete spatial randomness relies on comparing an 
observed spatial pattern with patterns expected for a 
homogeneous Poisson process. Measures of dispersion 
may be calculated from the distribution of individuals 
among sample plots or distances between individuals 
within a population (nearest neighbor distances; Diggle, 
1983). For example, we may calculate nearest neighbor 
distances for initial nesting attempts on a study plot, then 
simulate a random distribution with the same number of 
points within that area, and compare the observed nearest 
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with those 


completely random process (i.e., homogeneous Poisson). 


neighbor distances expected from a 
We plot observed nearest neighbor distances (empirical 
distribution function, solid line) against hypothetical 
nearest neighbor distances under complete spatial 
randomness (cumulative distribution function, straight 
dashed line) to evaluate whether an observed spatial 
pattern may have been generated by a completely 
random process (Fig 2; Dingle, 1983). If a spatial 
distribution is random there is a straight-line relationship 
between observed nearest neighbor distances and a 
completely random distribution (Fig. 2). This straight- 
line relationship is the probability that empirical nearest 
neighbor distances take on a value equal to or less than 
the nearest neighbor distances expected under complete 
spatial randomness. If our plot of observed nearest 
neighbor distances falls to the left of the straight line, this 
indicates there is a greater probability of having nearest 
neighbor distances less than that under complete spatial 
randomness, a spatially aggregated distribution (Fig. 2). 
In contrast, if our plot of empirical data falls to the right 
of the straight line, this indicates fewer short distances 
than expected under complete spatial randomness, a 
spatially even distribution (Fig. 2). 

We used the SAS macro collection for nearest 
neighbor analysis of a spatial point pattern (NNASPP, 
Schabenberger & Pierce. 2002) to determine whether 
Western Sandpiper initial nesting attempts were spatially 
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Fig. 2 Three hypothetical spatial distributions (random, aggregated, and even) with corresponding plots of observed nearest 
neighbor distances along the y-axes (empirical distribution functions, solid line) and nearest neighbor distances expected 
under complete spatial randomness along the x-axes (cumulative distribution function, straight dashed line). Irregular 
dashed lines in each plot provide upper and lower simulation envelopes for nearest neighbor distances under complete 
spatial randomness and represent sample averages. 
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aggregated on our study plot annually. We modified two 
macros within this collection to account for the shape of 
our study plot. We altered the %unitco macro so that it 
did not scale nest coordinate data to a unit square; 
instead, this macro plotted nest locations within the 
specified domain of our 36 ha study plot. We also 
updated the %hpp macro to simulate a homogenous 
Poisson processes within the domain of our 36 ha study 
plot. This macro collection computed nearest neighbor 
distances for initial nesting attempts on the 36 ha plot 
each year, simulated 1000 datasets for each year with the 
same number of nests randomly distributed across the 
plot to provide an estimate of complete spatial 
randomness, calculated sample averages for estimates of 
complete spatial randomness (upper and lower 
simulation envelopes), and compared observed nearest 
neighbor distances with estimates of complete spatial 
randomness using Monte-Carlo significance tests. We 
also used the SAS macro collection NNASPP to 
determine if displaying males were spatially aggregated 
on our study plot on a daily basis during 2003-2005. 
gathered by two 


simultaneously surveying from the middle of the plot 


These data were observers 
towards its periphery. Care was given to avoid double 
counting of individuals while surveying the plot. We 
only recorded the initial location where displaying males 
were observed. For spatial analysis of displaying males, 
we selected all days between 10 May and 10 June for 
which we observed at least 20 displaying males on the 














plot “average number of days per year + SD = (17 + 2), 














mean number of displaying males per day + SD = (27 + 
3)” . We did not observe 220 males displaying on the 
plot in a single day outside of this period. 
1.4 Nearest neighbor distance and individual male 

behavior 

We examined individual male behavior and their 
spatial distribution prior to being paired with a fertile 
female during 2003—2005. We identified males that were 
not yet paired with fertile females using nest initiation 
data (Lanctot et al, 2000). From those, we selected males 
with a minimum of 10 display flight observations for 
analyses (minimum of four hours between consecutive 
observations). This process resulted in the selection of 32 
males over a three-year period with a mean of 18 
observations/individual. Using these data, we estimated 
relative display area size for each male by drawing 
minimum convex polygons around the location of 
display flights. We also computed the proportion of 
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observations in which males exhibited behavior in five 
behavioral classes (display flight, chase, fight, solicit, 
and self-maintenance). We considered a male to have 
performed a display flight if he vocalized at a height of at 
least 5 m for a minimum of 10 s, chasing behavior 
entailed a male flying after one or more conspecifics, 
physical contact with a conspecific denoted fighting 
behavior, the tail-up courtship stance described by 
Holmes (1973) was used to quantify solicitation, and 
self-maintenance behaviors included loafing, sleeping, 
and preening. As a measure of spatial dispersion, we 
calculated the distance from each initial nesting attempt 
on the plot to its two nearest neighbors using ArcMAP. 
We then examined the relationship between male 
behaviors prior to being paired with a fertile female and 
the sum of the distances from each male’s initial nesting 
attempt to the two nearest neighboring nests using 
Spearman rank correlation (Zar, 1999). 
1.5 Testing the Material Resources hypothesis 
1.5.1 Habitat use in relation to availability 

To examine the relationship between spatial 
dispersion and habitat patches (tundra, tundra-grass, 
undulating-tundra), we delineated aggregated and 
dispersed nesting birds by separating nearest neighbor 
distances (sum of the distance from each nest to its first 
two nearest neighbors) into quartiles annually. This 
resulted in nearly a three-fold difference in nearest 
neighbor distances between the first and fourth quartiles 
or aggregated and dispersed nesting birds (Tab. 1). We 
then overlaid initial nesting attempts on our vegetation 
map (Fig. 1), and calculated the number of nests within 
each habitat class to determine whether aggregated (1“ 
nearest neighbor distance quartile) and dispersed (4" 
nearest neighbor distance quartile) nesting birds 
differentially nested in specific habitat patches. We 
compared the number of aggregated and dispersed 
nesting birds in each habitat type annually to the number 
of nests expected in each habitat type if birds utilized 
habitat in relation to availability using chi-square 
goodness of fit tests (Zar, 1999). Similarly, we evaluated 
the locations where the earliest arriving males (first five 
days males were observed displaying each year) were 
observed performing display flights. We compared the 
number of displaying males in each habitat patch type to 
the number of males expected in each habitat type if 
birds utilized habitat in relation to availability using chi- 
square goodness of fit tests. We also examined male and 
female age in relation to initial nest placement among the 
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Tab. 1 Sum of distances between initial Western Sandpiper 
nesting attempts and two nearest neighbor nests on 
a 36ha plot at Kanaryarmiut Field Station, Yukon- 
Kuskokwim Delta, Alaska (1999-2001, 2003-2005). 











Quartile (N) Median (m) Mean (m) + SD Range (m) 
1 (120) 52 50+15 13-74 
2 (120) 75 7549 53-97 
3 (117) 99 98 +12 75-124 
4 (121) 137 147 + 36 102-278 
three habitat patch types (tundra, tundra-grass, 


undulating-tundra). We grouped birds into three age 
classes for this analysis (minimum age, number of years 
observed breeding at the site = 1, 2, 23). We randomly 
selected a single nesting observation for birds observed in 
multiple years for this analysis to avoid pseudoreplication, 
and used chi-square goodness of fit tests to determine if 
habitat patch use among age classes varied from that 
expected based on availability. 
1.5.2 Microhabitat associated with nests 

During 2004, we examined microhabitat associated 
with Western Sandpiper nests (N = 128) by comparing 
the proportion of five vegetation types (birch, grass, herb, 
moss, lichen) and the amount of vertical relief within a 1 
m° frame centered on each nest and a random point 
within 40 m of each nest. We used matched-pairs logistic 
regression to examine microhabitat variation between 
nest sites and associated random points (Johnson & 
Oring, 2002). We compared potential models beginning 
with a model containing all six predictor variables and 
compared this model’s fit with simpler models to achieve 
the most parsimonious regression model that accurately 
represented the data. Akaike information criteria (AIC) 
were computed for each model to compare model ranks 
and weights in determining the likelihood of each model 
given the data (Burnham & Anderson, 1998; Anderson et 
al, 2000). 
1.6 Testing the social mate choice hypothesis 
1.6.1 Minimum age in relation to nest dispersion 

We used contingency tables to examine the 
relationship between Western Sandpiper minimum age 
and nest dispersion (Zar, 1999). We constructed a 2 x 3 
contingency table for each sex with aggregated (1“ 
nearest neighbor distance quartile) and dispersed (4" 
nearest neighbor distance quartile) nesting birds as rows 
and minimum age (number of years observed breeding at 
the site = 1, 2, 23) as columns to test whether Western 
Sandpiper age was independent of nest dispersion. 


1.6.2 Clutch initiation date in relation to nest dispersion 
and habitat 

We used SAS PROC MIXED to compare mean 
clutch initiation date, for initial nesting attempts between 
aggregated and dispersed nesting birds using analysis of 
variance (ANOVA, Littell et al, 2002). Initiation date 
was the dependent variable in the model, and nest 
dispersion was a two-level fixed factor independent 
variable (aggregated = 1“ 
quartile, dispersed = 4" 


nearest neighbor distance 
nearest neighbor distance 
quartile). We also included year as a random factor in the 
model to control for annual variation. Residuals were 
normally distributed (Shapiro-Wilk tests P > 0.10) and 
had homogeneous variance after a square root 
transformation. We compared mean initiation date using 
Least Square Means (a = 0.05). Based on habitat 
analyses, we also performed an a posteriori analysis 
comparing mean clutch initiation date between birds that 
nested in or out of undulating-tundra habitat while 
controlling for year effects. This analysis was similar to 
that described above except the independent variable was 
a two-level fixed factor (nests in undulating-tundra 
habitat, nests not in undulating-tundra habitat). 
1.7 Testing the predation hypothesis 

We compared nest success between aggregated and 
dispersed nesting birds using logistic regression. For this 
model, the dependent variable, nest success (at least one 
egg hatched in a clutch), was a binomial response (1 = 
successful, 0 = unsuccessful). Nest dispersion was a 
categorical variable with two levels (aggregated = 1“ 
nearest neighbor distance quartile, dispersed = 4" nearest 
neighbor distance quartile), and we included year 
(categorical) and initiation date (continuous) as 
covariates in the model as these factors are known to 
influence nest success in many avian species (Lack, 1954; 
Martin, 1987). Based on habitat analyses, we also 
performed an a posteriori analysis examining nest 
success between habitat types. As described above, we 
used logistic regression to model the probability of nest 
success based on whether birds nested in or out of 
undulating-tundra habitat while controlling for initiation 
date and year effects. In this model, habitat type was a 
categorical variable with two levels (undulating-tundra, 
not undulating-tundra). 


2 Results 


During six years of study (1999-2001 and 
2003-2005), we individually color banded 453 adult 
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Western Sandpipers and located and monitored 680 nests. 
Four hundred and thirty-three of these nests were initial 
nesting attempts on the 36 ha plot (mean number of 














initial nesting attempts per year + SD = (72 + 10). 
2.1 Aggregation patterns 

Displaying males exhibited an aggregated spatial 
distribution, and subsequent initial nesting attempts 
attended by those males also were aggregated. Displaying 
males exhibited an aggregated spatial distribution on the 
study plot during >80% of observation days when we 
observed = 20 displaying males (2003, 83% of 
observation days, NV = 18d, MC significance tests P < 
0.01 = 10d, P < 0.05 = 5d, P > 0.05 = 3d, 1000 
simulations/d; 2004, 82% of observation days, N = 17d, 
MC significance tests P < 0.01 = 11d, P< 0.05 = 4d, P > 
0.05 = 3d, 1000 simulations/d; 2005, 87% of observation 
days, N = 15 d, MC significance tests P < 0.01 = 10d, P 
< 0.05 = 3d, P > 0.05 = 2d, 1000 simulations/d). Western 
Sandpiper initial nesting attempts were spatially 
aggregated on the 36ha plot throughout the course of this 
study (1999, 7 = 31.0, P < 0.001, N = 49; 2000, 7 = 
36.2, P< 0.001, N = 84; 2001, 7 = 35.2, P< 0.001, N= 
85; 2003, 7 = 41.8, P < 0.001, N = 69; 2004, 7 = 38.0, 
P < 0.001, N = 83; 2005, 7 = 39.0, P < 0.001, N = 63; 
1000 simulations/yr). The observed distribution of 
nearest neighbor distances for initial nesting attempts in 
2003 (Fig. 1) indicated a spatially aggregated distribution 
(Fig. 3). The plot of observed nearest neighbor distances 
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Fig. 3 Plot of the empirical distribution function of observed 
nearest neighbor distances for Western Sandpiper 
initial nesting attempts on the study plot in 2003 (solid 
line) versus the cumulative distribution function of a 
hypothetical distribution of nearest neighbor distances 
under complete spatial randomness (straight dashed 
line). Irregular dashed lines provide upper and lower 
simulation envelopes for the cumulative distribution 
function based on 1000 Monte Carlo simulations. 
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(solid line, Fig. 3) falls to the left of the straight line 
relationship expected under complete spatial randomness 
(straight dashed line, Fig. 3) and is outside the simulation 
envelope for complete spatial randomness (irregular 
dashed lines, Fig. 3) indicating significant spatial 
aggregation. For brevity, we present the map of initial 
nesting attempts and plot of nearest neighbor distribution 
against complete spatial randomness for a single year 
(2003; Fig. 1,3); the comparable results for the remaining 
five years of study were similar. 
2.2 Nearest neighbor distance and individual male 

behavior 

Dispersed nesting males performed song-flight 
displays over a relatively larger area compared to 
aggregated nesting males prior to being paired with a 
fertile female (Fig. 4). Dispersed nesting males also were 
observed chasing conspecifics and soliciting females 
more often compared to aggregated nesting males (Fig. 
4). There was no correlation between nearest neighbor 
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Fig. 4 Spearman rank correlations between individual male (V 
= 32) behaviors (proportion of observations observed 
soliciting females and chasing conspecifics, and display 
area size) prior to being paired with a fertile female and 
the sum of the distances from subsequent initial nesting 
attempts to the two nearest neighbor nests (2003—2005). 
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distances and the number of song-flight displays, amount 
of fighting, or self-maintenance behavior observed (all r 
< 0.2, all P > 0.25, N= 32). 
2.3 Material resources hypothesis 
2.3.1 Habitat use in relation to availability 

During 2003-2005, 
displayed more often over tundra habitat and less often 


the earliest arriving males 


over undulating-tundra habitat than expected if birds 
utilized habitat patches in relation to availability (Fig. 5). 
Among initial nesting attempts, dispersed nests were 
more often in tundra habitat patches, less often in tundra- 
grass, and occurred in relative proportion to availability 
of undulating-tundra (Fig. 6). Whereas, aggregated initial 
nesting attempts were observed more often in 
undulating-tundra habitat patches, less often in tundra, 
and in relative proportion to availability in tundra-grass 
(Fig. 6). Younger males were more often associated with 
nests in undulating-tundra habitat patches, less often in 
tundra, and in relative proportion to availability in 
tundra-grass (Fig. 7). In contrast, older males were more 
often associated with nests in tundra habitat patches, less 
often in tundra-grass, and were never observed nesting in 
undulating-tundra (Fig. 7). Regardless of age, initial 
nesting attempts in the three habitat patches did not vary 
from that expected if females utilized habitat patches in 
relation to availability (lyr Loasa = 2.1, P > 0.25, N= 
89, 2yr Yo0s,2 = 1.1, P > 0.50, N= 27, S3yrs Z oos,2 = 
0.6, P > 0.50, N = 13). 
2.3.2 Microhabitat associated with nests 
Matched-pairs logistic regression revealed the 
proportion of birch (all models, 77; =11.1, P < 0.001) 
and graminoids (four top-performing models, 7°; =7.4, 


P < 0.01) within 1 m° (all other predictor variables, 7°, 
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Fig. 5 Proportion of male Western Sandpiper display-flights 
observed over three habitat patch types during the first 
five days males were observed displaying on the study 
plot annually (2003-2005) 

Habitat use by displaying males differed significantly from that expected if 
birds utilized habitat patches in relation to availability (2003 Zoss, 2 = 5.4, 
P < 0.10, N= 29; 2004 05,2 = 6.0, P < 0.05, N = 26; 2005 77o05,2 = 6.4, 
P < 0.05, N= 45). 
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Fig. 6 Proportion of initial nesting attempts by aggregated and 


% of nests in each habitat patch type 








dispersed nesting birds in three habitat patch types on a 

36 ha plot annually (1999-2001, 2003—2005) 
Habitat use for both aggregated and dispersed nesting birds differed 
significantly from that expected if birds utilized habitat patches in relation to 
availability (* = P < 0.05). Aggregated nesting birds: 1999 700s, = 6.3, P 
< 0.05, N= 12, 2000 005,2 = 6.2, P< 0.05, N= 23, 2001 7o05,2 = 14.3, P 
< 0.001, N= 26, 2003 77005,2 = 8.2, P< 0.025, N= 17, 2004 7o0s,2 = 3.0, 
P> 0.10, N = 23, 2005 7 005,2 = 6.3, P < 0.05, N = 16; Dispersed nesting 
birds: 1999 7o0s,2 = 2.3, P > 0.25, N=17, 2000 705.2 =6.6, P< 0.05, N 
= 22, 2001 7o05,2 = 7.0, P < 0.05, N= 24, 2003 7’o05,2 =5.1, P< 0.10, N 
= 17, 2004 7o05,2 = 6.5, P < 0.05, N= 20, 2005 7’o05,2 =5.9, P< 0.05, N 
=18. 
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Fig. 7 Proportion of initial nesting attempts in three habitat 
patch types in relation to minimum male age (# yrs 
observed breeding at the site) 

Habitat use for the youngest and oldest age classes (1* year and > 3 years, 
respectively) differed significantly from that expected if birds utilized 
habitat patches in relation to availability (* = P < 0.05, lyr Kos, 2 =12.2,P 
< 0.01, N=201, 2yr 7’o0s,2 = 0.6, P > 0.95, N= 20, > 3yrs 7o05,2 = 6.8, P 
< 0.05, N= 17). 


<1.0, P > 0.5, Tab. 2) were the dominant microhabitat 
predictors of nest presence. The top-performing model 
contained only birch and graminoids and this model had 
> 2 times the support of the next best performing 
candidate model (Tab. 2). There was a consistent 
positive association of birch and a negative association 
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Tab. 2 Multi-model inference based on results from matched-pair logistic regression of Western 
Sandpiper nest sites and matched random points 
Odds ratio 
Model AIC Ai Wi 
Birch Graminoids 
birch, graminoids 154.62 0.0 0.58 1.08 0.95 
birch, graminoids, herbs 156.22 1.6 0.26 1.08 0.95 
birch, graminoids, herbs, relief 158.12 3.5 0.10 1.08 0.95 
birch, graminoids, herbs, relief, moss, lichen 160.02 5.4 0.04 1.09 0.96 
birch 161.79 7.2 0.02 1.08 — 


Model parameters denote the proportion of each vegetation type within 1m? of nest sites and matched random points. 4; is the rank of each model by rescaling 


AIC values such that the model with the minimum AIC value has a value of zero, and w; is the likelihood of the model given the data. Odds ratios indicate 


percentage increase (birch) or decrease (graminoids) in the odds of a nest being present (i.e. 8%—9% increase in the odds of a nest being present when birch was 


present, 4%—5% decrease in the odds of a nest being present when graminoids were present). 


of graminoids with nest sites compared to random points 
in the four top-performing models (Tab. 2; range of odds 
ratios, birch = 1.08—1.09, graminoids = 0.95—0.96). 
2.4 Social mate choice hypothesis 
2.4.1 Minimum age in relation to nest dispersion 

Results did not indicate a difference in male or 
female age between aggregated and dispersed nests 
(males 7os,2 = 3.10, P> 0.10, N = 238, females 70s, 2 
= 1.59, P > 0.25, N= 129). 
2.4.2 Clutch initiation date in relation to nest dispersion 

and habitat 

Mean clutch initiation date was June 1* “mean 
Julian date + SE all years combined = (152 + 2)” . 

There was no difference in mean clutch initiation 








dates between aggregated and dispersed nests (f))3 = 
—1.50, P = 0.14), or between nests located in or out of 
undulating-tundra habitat (t425=—0.91, P = 0.36). 
2.5 Predation hypothesis 

Logistic regression did not indicate that nest success 
significantly varied between aggregated and dispersed 
nesting birds (7, = 0.02, P = 0.90, N = 220); however, 
there was significant annual variation in nest success 
during the course of study (^s = 11.76, P = 0.04, N = 
220). The annual percentage of aggregated nests that 
hatched ranged between 1%-—36% and the percentage of 
dispersed nests that hatched annually ranged between 
11%—43%. The logistic regression model examining nest 
success between habitat types indicated that nest success 
varied between nests located in and out of undulating- 
tundra habitat (Fig. 8, 77, = 4.2, P = 0.04, N = 433). 
Nests that were not in undulating-tundra patches were 
63% more likely to hatch compared to nests within 
undulating-tundra patches (odds ratio = 1.63, 95% Wald 
confidence limits 1.02-2.61). 
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Fig. 8 Percentage of Western Sandpiper initial nesting attempts 





to hatch (+ SD) in two habitat categories annually on a 
36ha study plot at Kanaryarmiut Field Station, Yukon- 
Kuskokwim River Delta, Alaska 


Numbers above bars indicate sample sizes. 


3 Discussion 


Western Sandpipers exhibited aggregated breeding 
behavior at a relatively small spatial scale (36 ha). Prior 
to clutch initiation, male song flight displays were 
spatially aggregated and subsequent initial nesting 
attempts on the plot also were. Although many taxa of 
non-colonial birds have been reported to aggregate 
breeding territories (Stamps, 1988; Reed & Dobson, 
1993; Cade & Woods, 1997; Danchin & Wagner, 1997; 
Reed, 1999; Etterson, 2003;Tarof et al, 2004), to our 
knowledge, this is the first study to quantify aggregated 
breeding behavior in a non-colonial monogamous 
shorebird species prior to hatch. 

3.1 Social mate choice hypothesis 

If sexual selection favors male aggregation because 
it facilitates social mate choice (Allee, 1951; Darling, 
1952), then breeding aggregations of Western Sandpipers 
may represent an adaptive male strategy. Facilitation of 
mate choice would be especially important among birds 


No. 2 Matthew Johnson et al: Proximate and ultimate factors that promote aggregated breeding in the Western Sandpiper 137 


in which female choice is constrained temporally 
because of a compressed breeding season (Veen et al, 
2001) or when male parental care is important (Gowaty, 
1996), such as in the Western Sandpiper. 

We would expect aggregated males to acquire a 
mate (and subsequent clutch) earlier than more dispersed 
individuals if aggregated males are more attractive than 
solitary males, or if unpaired females more easily detect 
aggregated males (Real, 1990). We failed to find support 
for the social mate choice hypothesis in Western 
Sandpipers as mean clutch initiation date did not vary 
between aggregated and dispersed nesting birds. Because 
of the time constraints associated with arctic breeding, it 
is possible that there is little-to-no opportunity for 
substantial variation in clutch initiation date for Western 
Sandpipers (Johnson & Walters, 2008; Johnson et al, 
2010). Regardless of potential temporal constraints on 
clutch initiation date, male aggregations could still 
facilitate social mate choice if aggregated males were 
able to attract higher quality females compared to 
dispersed individuals. However, we found no evidence 
that more experienced females were associated with 
aggregated nests compared to dispersed nests. 

3.2 Predation hypothesis 

We failed to find support for the predation 
hypothesis as an explanation for aggregated breeding in 
Western Sandpipers. The predation hypothesis (Hamilton, 
1971; Pulliam, 1973) predicts that bird’s aggregate 
because clustering of individuals reduces predation via 
proximate antipredator mechanisms. However, results 
from six years of data revealed no consistent variation in 
nest success between aggregated and dispersed nesting 
birds. 

3.3 Material resources hypothesis 

The material resources hypothesis predicts that 
individuals aggregate in response to patchily distributed 
resources (vegetation and/or food, Kiester & Slatkin, 
1974). We found support for the material resources 
hypothesis in that aggregated nests were more often 
associated with undulating-tundra habitat patches and 
dispersed nests were more often in tundra patches. The 
disproportionate use of undulating-tundra patches by 
aggregated birds resulted in a reduction in their use of 
tundra patches, but they did place nests in tundra-grass 
patches in relative proportion to habitat availability. In 
contrast, as a result of their disproportionate use of 
tundra patches, dispersed nesting birds used tundra-grass 
habitat patches less than expected, but used undulating- 


tundra habitat 
availability. 


patches in relative proportion to 

Previous investigations have suggested that social 
organization, territorial behavioral, and nest dispersion 
have coevolved to increase a species’ ability to exploit 
resources (food, habitat, and mates) over both space and 
time (Holmes, 1966, 1972; Oring & Knudson, 1972; 
Graul, 1973; Pitelka et al, 1974; Safriel, 1975; Emlen & 
Oring; 1977). Among monogamous Arctic breeding 
sandpipers, such as the Western Sandpiper, it has been 
assumed “...that most if not all species in this group are 
dispersed relatively evenly over the available habitat 
each year.” (Pitelka et al, 1974). Although territory size, 
and subsequently breeding density, may vary as a result 
of whether food is primarily acquired on (Dunlin, C. 
alpina) or off (Western Sandpiper) the nesting territory 
(MacLean, 1969; Soikkeli, 1967; Holmes, 1970, 1971), 
large fluctuations in densities within a species at 
particular sites have not been reported (Pitelka et al, 
1974). Holmes (1971) suggested that high densities of 
Western Sandpipers in some areas of western Alaska are 
related to the patchy distribution of suitable nesting 
habitat (upland tundra; dwarf shrub-heath tundra in 
Holmes (1971)). As a result, Holmes (1971) reported that 
the overall distribution of Western Sandpiper in western 
Alaska is clumped, but within suitable nesting habitat 
birds are regularly dispersed. Our results are contrary to 
this observation. Although Holmes’ study site was 
merely 30 km northwest of ours, it is possible that upland 
tundra vegetation was more homogeneous at his site 
compared to ours, and this could result in the 
discrepancy in observed breeding distribution between 
studies. Holmes (1971) did not delineate vegetation 
patches within upland tundra habitat as we did, thus we 
are not able to compare upland tundra vegetation 
community composition between the two sites. 
3.4 Ideal free and ideal despotic distributions 

Our results can be interpreted in terms of the theory 
of ideal free and ideal despotic distributions (Fretwell & 
Lucas, 1970; Fretwell, 1972). In situations conforming to 
the ideal free distribution, animals move freely among 
habitats and assort themselves in proportion to resource 
availability. However, once a critical density is attained 
in preferred habitats, individual fitness is reduced in 
preferred habitats and individuals begin to colonize less 
preferred habitat where competition is less. This results 
in individual fitness being equal over a range of habitats, 
resources, or other conditions. Alternately, the ideal 
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despotic distribution model predicts that subordinate 
individuals are constrained in their choice of habitat by 
dominant individuals (Messier et al, 1990), resulting in 
differences in fitness among habitats. Primary 
assumptions of this model are that different habitats 
provide varying opportunities for individual fitness, and 
that competition will ensue where resources are limited 
(Leibold, 


dispersed individuals more often nested in tundra habitat 


1995). Based on our observations that 


patches (Fig. 6), displayed over larger areas, chased 
conspecifics and solicited females more often (Fig. 4), 
that the earliest arriving males primarily displayed over 
tundra and tundra-grass habitat patches (Fig. 5), and that 
older males nested in tundra and tundra-grass but not in 
undulating-tundra patches (Fig. 7), we conclude that 
there is competition among males for tundra and tundra- 
grass habitat patches. 

For the breeding distribution of Western Sandpipers 
to conform to an ideal deposit distribution, fitness must 
vary among habitat patches. On our study plot, nests in 
tundra and tundra-grass habitat patches that males 
appeared to prefer were more likely to hatch compared to 
nests in undulating-tundra patches (Fig. 8). We can only 
speculated as to why nest success was lower in 
undulating-tundra patches, but we did observe mink 
(Mustela vison) regularly using these patches to move 
across the landscape from one lakeshore to another, and 
both Arctic (Alopex lagopus) and red (Vulpes vulpes) fox 
commonly used vertical structure, such as frost heaves, 
to apparently surprise prey while foraging near the study 
site. Greater vertical relief associated with undulating- 
tundra patches also may provide such concealment for 
foraging fox on our study plot. Regardless of the 
mechanism(s) underlying variation in reproduction, we 
conclude that Western Sandpipers exhibited a despotic 
breeding distribution on the study-plot, with subordinate 
males forced to aggregate their nests in undulating- 
tundra habitat patches where fitness is reduced. We 
typically think of breeding aggregations as being 
concentrated on or near a critical resource, Western 
Sandpiper breeding aggregations appear to result from 
dominant and/or older individuals excluding younger, 
subordinate individuals from preferred habitat. 

3.5 Microhabitat associated with nests 

There is one intriguing inconsistency in our 
interpretation of our results. Dispersed males nested in 
tundra habitat patches more often than expected based on 
availability, less often in tundra-grass than expected, but 
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in relative proportion to availability of undulating-tundra 
patches. If nesting in undulating-tundra resulted in 
reduced fitness via lower nest success, why were older, 
more aggressive, males utilizing undulating-tundra at all? 
Why were these individuals not utilizing tundra and 
tundra-grass patches to the exclusion of undulating- 
tundra? Analysis of nest site microhabitat features 
revealed that Western Sandpiper nests were more often 
associated with areas containing a higher percentage of 
dwarf birch and a lower percentage of graminoid species. 
It is reasonable to suppose that reduced use of tundra- 
grass patches by dispersed nesting birds, but continued 
use of undulating-tundra patches, is the result of birds 
avoiding graminoid species and preferentially placing 
nests near dwarf birch. However, such reasoning does 
not explain nest placement among aggregated nesting 
birds, as those birds used tundra-grass habitat patches in 
relative proportion to availability, but reduced their use 
of tundra patches. Continued use of tundra-grass and 
reduced use of tundra patches by aggregated individuals 
may be the result of aggregated individuals making the 
best of a bad situation. Aggregated males appeared to be 
subordinate to dispersed individuals, so they may have 
been excluded from using habitat other than tundra-grass 
and undulating-tundra patches. 


4 Conclusions 


Western Sandpipers exhibited a spatially aggregated 
breeding distribution on a 36 ha plot. Displaying males 
were aggregated on a daily basis as were initial nesting 
attempts on an annual basis. We found support for the 
material resources hypothesis, as dispersed individuals 
were more often associate with tundra habitat patches, 
and aggregated individuals nested more often in 
undulating-tundra patches. There also was support for an 
ideal despotic distribution with lower reproductive 
success associated with aggregated nesting birds in 
undulating-tundra patches. Although we may typically 
think of breeding aggregations as being concentrated on 
or near a critical resource, Western Sandpiper breeding 
aggregations appear to result from dominant and/or older 
individuals excluding younger, subordinate individuals 
from preferred habitat. 
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